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Abstract 

We show that the transverse momentum dependent transversity function is 
proportional to the longitudinal polarization of a quark in a transversely polarized 
proton. This result suggests an alternative, convenient method for determining 
transversity, without knowing unusual fragmentation functions. The method 
consists of measuring the double spin azimuthal asymmetry in semi-inclusive pion 
leptoproduction by a transversely polarized proton target. The asymmetry, which is 
twist 3, is estimated to be more than 10% under the most favorable conditions. The 
experiment we suggest is feasible at facilities like DESY and CERN. 



PACS numbers: 13.85.Qk, 13.88.+e 



1 Introduction 



The transversity function, usually denoted as hi{x), is one of the three leading 
twist distributions of the quark inside the nucleon, together with the unpolarized 
distribution, q{x), and with the longitudinally polarized one, Aq{x). In an infinite 
momentum frame where the nucleon spin is perpendicular to its momentum, one has 



where q^^i) {x) is the probability density to find a quark whose spin is parallel (oppo- 
site) to the nucleon spin. Important information on quark dynamics inside the nucleon 
can be extracted from transversity. For example, the difference hi{x) — Aq{x) is quite 
sensitive to the orbital angular momentum of the quarks. This is why high energy 
spin physicists have been concentrating their efforts for several years in determin- 
ing transversity [1-5]. However, this appears a particularly difficult task[6-10], owing 
to the chiral odd character of hi, unlike q and Aq. Difi^erent observables sensitive 
to transversity have been singled out. Among them we recall the Drell-Yan (DY) 
double spin asymmetry [1, 9] and the single [7, 11, 12, 13] and double [14, 15] spin 
asymmetry in semi- inclusive deep inelastic scattering (SIDIS). 

For the moment SIDIS reactions - realized or planned by HERMES[11, 12] , SMC[13] 
and COMPASS collaborations - appear most promising for determining hi. These 
reactions are of the type 



where is a charged lepton, p a polarized proton target and h a hadron. If h 
is a spinning particle whose polarization can be determined (e. g., a A), reaction 
(2) is kinematically isomorphic to the double spin DY process. However, as we shall 
discuss in the following, it is more convenient to detect a pion in the final state. In 
this case we are faced with a single spin asymmetry, therefore we have to exploit a 
possible azimuthal asymmctry[4, 16, 17, 18, 19, 20] of reaction (2). Indeed, taking 
a frame whose origin coincides with the location of the proton target and whose z- 
axis is along the momentum transfer of the lepton, the angular distribution of the 
pion may present, in principle, an azimuthal asymmetry with respect to the plane 



hi{x) = q^{x) - qi{x), 



(1) 



Ip^fhX, 



(2) 
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passing through the proton spin and through the 2;-axis. Incidentally, such kind of 
asymmetry may be sensitive to various transverse momentum dependent (t.m.d.) 
distributions [16, 17, 21, 22, 23, 24], defined as "new" by Kotzinian and Mulders[21] 
(KM); see also Mulders and Tangerman[16] (MT). In order to extract distributions 
from azimuthal asymmetries, it is convenient to define weighted asymmetries [19, 21], 
that is, ratios of the type 

, , . f wd(j 

{K) = (3) 

Here w is a given function of the azimuthal angle of the final hadron h, over which 
we perform the integrations indicated in formula (3). In the reaction considered the 
azimuthal asymmetry is caused by the Collins effect [25], consisting of an interference 
term[26]. This effect is described by a T-odd (and chiral odd) t.m.d. fragmenta- 
tion function, say k. The corresponding weighted asymmetry results to be sensitive 
to the product hi{x)c{z), where c — J wndcf) is the so-called Collins fragmentation 
function [25]. As claimed by Jaffe[4], this may become the "classic" way of determin- 
ing the proton transversity distribution, provided c{z) is known to some precision and 
is not too small. But at present we know very little about this function, moreover the 
methods suggested for inferring it from data are complicated and require an adequate 
statistics[20, 4]. 

Analogous considerations could be done about the method proposed by Jaffe and 
Ji[27] (JJ). This consists of measuring the double spin asymmetry in a SIDIS reaction 
of the type (2), where the initial lepton is longitudinally polarized, the proton target 
is transversely polarized and the final hadron is a pion, i. e., 

ip^ f ttX. (4) 

The corresponding asymmetry is defined as 

^(|k|;ft.;n,) = ^^;^^^^;^ . (5) 

Here, as usual, u is the lepton energy transfer and = —q^, q being the four- 
momentum transfer. Furthermore k is the momentum of the initial lepton and Hy 
the component of the pion momentum along the momentum transfer. Lastly da^-^ and 



d(7-f<_ are the polarized differential cross sections for reaction (4) integrated over the 
azimuthal angle of the pion, arrows indicating the proton and lepton polarization. 
Asymmetry (5) is sensitive to the product hi{x)e{z)[27], where e{z) is the twist-3 
fragmentation function of the pion. The extraction of hi depends again critically on 
an unknown function. 

If the cross section is not integrated over the transverse momentum of the final 
pion, reaction (4) [21, 17] exhibits an azimuthal asymmetry. This is sensitive[21] 
to the "new" function qit, proportional to the longitudinal quark polarization in a 
transversely polarized proton. 

The aim of this paper is to re-examine such an azimuthal double spin asymmetry. 
We derive the differential cross section for reaction (2), starting from the definition 
of t.m.d. transversity function as given by Jaffe and Ji[2] (JJl), i. e., 

5q^{x,p^)^ ^ 2TqT{x,p^). (6) 

T=±l/2 

Here, analogously to formula (1), qT{x,p±) is the probability density to find, in a 
transversely polarized proton, a quark whose spin is parallel {T — 1/2) or opposite 
(T = -1/2) to the proton spin. That is, instead of the usual helicity representation, 
we consider a canonical one, such that the quantization axis is taken along the pro- 
ton polarization. The asymmetry we calculate turns out to coincide with the one by 
KM, provided we identify * gir with Sq±. We shall prove this identity for massless 
quarks, which amounts to saying that, owing to transverse momentum, a quark in a 
transversely polarized proton has a longitudinal polarization, related to transversity. 
Therefore Sq± - denoted as hir by MT - plays a major role in the azimuthal asym- 
metry of reaction (4), where the target is transversely polarized. Moreover, as we 
shall see, this distribution is somewhat relevant also in the case of a longitudinally 
polarized target. All this suggests an alternative, convenient method for determining 
the transversity. Indeed, as a consequence of our result, 5q± contributes not only to 
the chiral-odd component of the t.m.d. correlation matrix[16, 23, 28], but also to its 
chiral-even part. Therefore this distribution may be coupled to a chiral-even fragmen- 
tation function, which is generally easier to determine than a chiral-odd one. This 
* up to a normalization constant, see appendix and subsect. 2.4 
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result is not completely surprising, since we have shown in a previous paper[29] that 
DY from singly polarized proton-hadron collisions - which is kinematically isomorphic 
to reaction (4) - produces a muon polarization sensitive to 5q^. 

Drawbacks, common to any SIDIS reaction, have been pointed out rcccntly[30, 
26, 31, 32]. In fact, the asymmetries measured in these reactions are likely to derive 
non-negligible contributions from spurious effects, caused by interactions [30, 26] or 
correlations [32] between the active quark and the spectator partons, or by fragmen- 
tation from the target remnant diquark, or by hght cluster decays[31]. Such effects 
have to be taken into account by means of suitable fracture functions[26]. 

Sect. 2 is dedicated to the derivation of the formulae for the cross section and for 
the asymmetry we are interested in. Moreover we illustrate some consequences of the 
identity Qit = in the chiral limit, which we prove in the appendix. In sect. 3 we 
discuss how to infer the transversity from data of azimuthal double spin asymmetry, 
taking into account the spurious effects. In sect. 4 we compare our method with 
others. Lastly sect. 5 is devoted to a short conclusion. 



2 Cross section and azimuthal asymmetry 
2.1 Cross section 

We calculate the differential cross section for reaction (4) in the framework of a 
QCD- improved parton model[33]. For the moment we neglect the contribution of the 
fracture function, which will be discussed in sect. 3. In the laboratory frame the 
differential cross section reads, in one-photon exchange approximation, 

where M is the proton rest mass and L^i, {H^i,) the leptonic (hadronic) tensor. dT, 
the phase space element, reads 

dr ^ j^d^k' 6{k") 9{k'o) d^P 6{P^-ml) 9{Po). (8) 



Here k' and P are, respectively, the four-momenta of the final lepton and of the pion, 
whose rest mass is 171^^. The leptonic tensor is, in the massless approximation, 

V = ^Tr[|^(l + A,75)7Ml^'7.], (9) 

and k = k' + q (/c^ = 0) being respectively the helicity and the four-momentum of 
the initial lepton. Trace calculation yields 

L^iu = kj^kl + k'^ky - g^^yk ■ k' + iXiEanpuk"'^'^ ■ (10) 

As regards the hadronic tensor, the generahzed factorization theorem[34, 23, 35] in 
the covariant formalism[36] yields, at zero order in the QCD coupling constant, 

= \ jZe) I dV,^f{p'-P)hiM\S), (11) 
= Y:dip)Tr{p'^l,p'l^^)■ (12) 

L 

Here the factor 1/3 comes from color averaging in the elementary scattering process 
and / runs over the three light flavors {u,d,s) and antiflavors {u,d,s), Ci = —64 = 
2/3, 62 — €3 — —65 = —66 = -1/3. p and p' are respectively the four-momenta of 
the active parton before and after being struck by the virtual photon. S is the Pauli- 
Lubanski (PL) four- vector of the proton, qiip) is the probability density of finding a 
quark (or an antiquark) with four-momentum p, whose spin is parallel {L = 1/2) or 
antiparallel (L = — 1/2) to the proton spin. Analogously ^p-^{p') is the fragmentation 
function of a quark of four-momentum p' into a pion of four-momentum P. Moreover 

dr, = j^^d'p 6{p')9{po)d'p' 5{p')e{p',) 5\p' -p-q), (13) 

the active parton being taken on shell and massless. Lastly the p's are the spin density 
matrices of the initial and final active parton, i. e.[29], 

p^ = i^[l + 2L75(A + ^)] and P' = (14) 

Here 2Lr) is the transverse PL four-vector of the active parton, while A is the longi- 
tudinal component of the quark spin vector. Formulae (14) are consistent with the 
Politzer theorem[37] in the parton model approximation. These imply, together with 
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eq. (12), that rj does not contribute to hj^^. For later convenience we re- write this 
last tensor as 

hL = 7 WiP)^t^'^ + A5g^(p)a J . (15) 



Here 

Si,u = rr(^7^/7^), a^^ = Tr{-tbinni>'lu)- (16) 

Moreover {p) = I]l=±i/2 9l (p) is the unpolarized quark distribution and 5q^{p) = 
X)l=±i/2 2Ivg{(p). A is a Lorentz invariant, pseudoscalar quantity, such that |A| < 1. 
If we neglect the parton transverse momentum, the only way of constructing such a 
quantity with the available vectors is 

Here we have exploited the fact that z/ is a Lorentz scalar and that in the laboratory 
frame q = (z^, q) and S = (0, S), where S is the proton spin vector, = 1. Ay can 
be viewed as the helicity of the proton in a frame where the proton is moving along 
q. In this connection we observe that the spatial direction of the virtual photon 
can also be defined covariantly by means of the four-momenta of the photon and of 
the proton [19]. In order to take into account the transverse momentum, we have to 
adopt a frame where the proton momentum is large in comparison to M[38]. To this 
end we consider the Breit frame - coincident with the one adopted by Feynman[39] 
- where the virtual photon has a four-momentum q = (0, q^), with Iq^l = Q. In 
this frame the proton momentum is — ^qs, therefore the active parton carries a 
momentum = ~^^b + P±, where, as usual, x — Q'^/2Mi/ is the longitudinal 
fractional momentum and the transverse momentum with respect to q^. We 
decompose the proton spin vector S into a longitudinal and a transverse component, 
i. e., 

S = A|,-^ + S^, Sx-q = 0. (18) 

|q| 

But the average spin of the quark is independent of the quantization axis, therefore 
the decomposition (18) implies 

XSq^ix, p±) = X\\SqUx, pi) + X±Sq{{x, (19) 



Here 



Ax = ^ (20) 



moreover Sq^^^x, pj_) (denoted as g\L{x, pi) by MT) is the t.m.d. helicity distribution. 
It is interesting to write the relationships of 5q^^ and 6q{^ to the usual distributions: 

Ag/(a;) = j cPpJqf^{x,^l)- h{{x) = I d^p±6qi{x,p^). (21) 

Now we carry on the integration (11) over the time and longitudinal components of 
p, taking the 2;-axis opposite to q. We get, in the light cone formalism, 

H,. = E 4 / d'p^if^iz, Pl)hi^{x, pr, S), (22) 

where the tensor hj^^ reads, after insertion of eq. (19) into eq. (15), 

= ^ {q^{x, pi)s^,^. + [\\Sq{\ {x, p1) + X±Sq{{x, a^^] . (23) 

q^ and (p^ are related, respectively, to the usual unpolarized distribution and to the 
pion fragmentation function, i. e., 

qf(x) = / d'p^qf(x, Pi), Dfiz) = / d'P^^f{z, Pi). (24) 

Moreover z — (J]| + Po)/(2|p'|) is the longitudinal fractional momentum of the pion 
resulting from fragmentation of the struck parton, whose momentum is p'. We have 



defined Pq = ^ml + P^, = P • p7|p'| and P^ = P - i^|p7|p'|, P being the pion 
momentum in the laboratory frame. Denoting by IIj^ the transverse momentum of 
the pion with respect to the photon momentum, we get 

P± = n^--2P±. (25) 

Therefore, if we keep 11^ fixed, P^ depends on p^. 

We notice that, although we have chosen a particular frame, the result (22) is 
apparently covariant. 
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2.2 Azimuthal asymmetry 

In order to calculate the double spin azimuthal asymmetry ^(|k|; Q, u; P) for reaction 
(4) - defined analogously to (5), but keeping the pion momentum P fixed - we have to 
substitute the leptonic tensor (10) and the hadronic tensor (22) into the cross section 
(7), taking into account relations (23) and (20). The result is 

T is the depolarization of the virtual photon with respect to the parent lepton[21]. 
Moreover we have introduced the quantities 

- g^(x,z,nl) = |dV9^(.^,pi)v'^(^,Pi), (27) 
SQ^ = (5gj[(a;,z,ni) + nx-S(5gi(x,^,nx), (28) 

5g{(x,^,ni) = A|||dV<^?f(^,Pi)<^^(^,Pl), (29) 
5gi(a;,;2,n^)n^-s = |dVA±5gi(x,px)<^^(^,Pi). (30) 

The pseudoscalar character of (eq. (28)) follows from assuming a massless lepton: 
indeed, the expression we have deduced for the asymmetry (see the first eq. (26)) 
holds in any frame where the lepton mass is negligible and the lepton helicity has 
the same value as in the laboratory frame. Below we shall show that bQ^ is twist 3. 
From formula (28) we deduce that, in order to maximize the contribution of to 
our asymmetry, one has to take the vector parallel to or opposite to S, that is, to 
select pions whose momenta lie in the (q, S) plane. Furthermore eq. (28) imphes that 
bQ^ is especially sensitive to (^?{(x, p_l) if q - S = 0. Indeed, in this situation the first 
term of eq. (28) - and more generally the JJ asymmetry - vanishes. Therefore events 
such that the lepton scattering plane is perpendicular to the proton polarization are 
particularly relevant to our aims. 

Since the products kj^k'_ and k-k'^ are invariant under boosts along the 2;-axis, 
we calculate them in the laboratory frame, where 

fc± = 1^(1 ± cos/3), fc; = ± cosie + /5)]. (31) 



Here k' = k — q is the final lepton momentum. Moreover /3 and d are, respectively, 
the angle between k and q and between k and k': 

|q|cos/3 = |k| - \\l\cose. (32) 

Now we consider the scaling limit, i. e., oo, v — > oo, jlMv — >■ x. Since 

(5^ 2|k||k'|(l — cos^), Q tends to zero in that limit, as well as 

-f^. ^-if^. y-^, (33) 

Then the second eq. (26) and eq. (17) yield, respectively, 

^ = 1 ^^2 ' '^ll ^ ~ — -sinOcoscI)', (34) 

l+(l-|/)2 

where (f)' is the azimuthal angle between the (k, k') plane and the (k, S) plane. There- 
fore 6Qp eq. (29), is twist 3, as follows from the second eq. (34) and from the first 
eq. (33). But also the second term of eq. (28) is twist 3, as is immediate to check. 
Therefore our asymmetry is twist 3. 



2.3 Numerical estimates 

Now we calculate the order of magnitude of the asymmetry (26) under optimal con- 
ditions. To this end, first of all, according to the considerations of subsection 2.2, we 
take into account events such that the azimuthal angle 0' (see the second eq. (34)) 
is about 7r/2 and IIj, is parallel (or antiparallel) to S. Moreover, eqs. (28) to (30) 
and eqs. (34) suggest that y should be chosen as close as possible to 1. Assuming the 
transverse momentum distributions involved in asymmetry (26) to be of the Gaussian 
type''", with equal widths of ~ 0.85 GeV/c, and setting |n_L| ~ 1 GeV and Q = 2.5 
GeV, the asymmetry (26) results in A ~ 0.4i?, where R = h{{x)/q^{x) has been 
determined by HERMES[11], \R\ = (50 ± 30)%. 



2.4 Remarks 

At this point some remarks are in order, 
"^see eqs. (37) to (39) below 
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(i) Parity invariance and a rotation of the reference frame by tt around the proton 
spin imply 

Sq± {x, P±)^Sq{ {x, (35) 

the typical property of a T-even function[40]. This, in turn, implies that, if we in- 
tegrate the SIDIS differential cross section - and therefore the hadronic tensor (22) 
- over n^, the second term of eq. (28) vanishes. On the other hand, upon in- 
tegration, the first term of eq. (28) goes over into X\\Aq^ {x)D-f (z), corresponding 
to the "kinematic" twist-3 term of the numerator in the J J asymmetry [27] . The 
above mentioned numerator includes also a "dynamic" twist-3 term of the type 
X\\[^h({x)^ef (z) + gj^{x)Df {z)][27] (see introduction), where gj'{x) is the transverse 
spin distribution. This term does not appear in asymmetry (26), since we limit our- 
selves to the contributions of the QCD parton model. 

(ii) Eqs. (26) to (30) and the first eq. (34) hold true for any orientation of 
the proton spin. However, if S is not perpendicular to the lepton beam. Ay does 
not decrease with Q. Therefore if, e. g., the proton is polarized longitudinally, the 
asymmetry derives still a contribution from the t.m.d. tranversity function, which, 
however, risks to be taken over by the JJ term[22]. We shall see in the next section 
a method for extracting 6qj_ under such unfavourable conditions, as, e. g., in the 
HERMES experiment [11]. 

(iii) The hadronic tensor (22), which we have derived starting from the definition 
of transversity by JJl, turns out to coincide with the tensor found by KM, provided 
we assume-!- 

9(t - '-^Sqi. (36) 

This relationship is proven in appendix, together with other useful connections among 
the t.m.d. distributions defined by MT (see also ref.[41]). Therefore KM's consider- 
ations and deductions about g(rp - like, e. g., its relation with g2 - can be applied to 
5qi. 

* according to the normalization constant chosen by KM for g(j. 
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3 Extracting transversity from data 

In this section we discuss how to extract Sq{^ from data. To this end we illustrate two 
different methods: a best fit to the asymmetry, with a suitable parametrization for 
Sq-j^, and the use of weighted asymmetries [21, 19]. Furthermore we suggest how to 
take into account the spurious effects, that is, contributions to asymmetry (26) not 
related to Sqj_. 

3.1 Gaussian parametrization 

A frequently used parametrization of the t.m.d. unpolarized distributions and frag- 
mentation functions consists of[16, 23, 35] 

q^{x,p]_) = {a/Tr)q^{x)exp{-ap'i), (37) 
ipfiz,Pl) = iajTr)Df{z)exp{-a^Pl). (38) 

Here a ~ 1.38 {GeV/c)~'^, as results from DY[42], and may be determined from 
two-jet events in e'^e" — > ttX. As regards the t.m.d. transversity, it looks appropriate 
to set [43] 

5q{(x, p_l) = {a/7r)h({x)exp{-ap\). (39) 
Here /i{ may be parametrized according to the suggestion of ref.[44], i. e., 

h({x) = l^^^^^^^ia + (3r^^ [qf{x) + Ag/(a;)] , (40) 
where N, a, (3 and b are free parameters, with |A^| < 1. 

3.2 Weighted asymmetries 
3.2.1 The KM weight function 

For the reaction we are studying, KM have defined a weighted asymmetry of the type 
(3), suggesting the weight function w{Il±) = 2Il± ■ S/M. Inserting this function and 
eq. (7) into eq. (3), and taking into account eqs. (10) and (22), we get 

E%ie}h{^'\x)Dfiz) 
E/=i ejqf{x)Df{z) 

12 



where 

h{^'\x) = ^ J d'p^ip^ ■ SfSql (42) 

The expression KM obtain for the weighted asymmetry differs from eq. (41) by the 
substitution hi {x) by QIt {x), where 

9(t\^) - ^ J d'p^pWiT- (43) 

f(l) f(l) 

But eq. (36) imphes gir (x) = hi '{x), so that the KM asymmetry turns out to 
coincide with eq. (41). Moreover the numerical evaluation of qit that KM infer from 
E143 data[45] allows, by means of eqs. (42) and (39), to evaluate hi. This behaves 
quite similarly to the bag model prediction (see, e. g., JJl), especially at small x. 
The accord between the two calculations can be made quantitative for Q of order 12 
to 15 GeV. 

3.2.2 Hcirmonic oscillator weight functions 

More refined information on the t.m.d. transversity function could be extracted from 
data by using a different set of weight functions, inspired to the Gaussian parametriza- 
tion. We assume for (p^ and Sqj_ expansions of the type 

<^^(^,Pi) =E<^k{z)^k{P): 5qi{x,p^) =J2ci{x)^K{p). (44) 

K K 

Here P = (P']_y^'^, p = (pi)^^^ and the $;^(P) are the eigenfunctions of the equation 



^k{P) = 0, (45) 



where ^ = id/dP. The fulfil an analogous equation, with a instead of a^^. These 
functions have been chosen in such a way that the terms with K — reproduce the 
parametrizations of subsection 3.1. Cj^{z) and c;^(x) are real coefficients. The set of 
weight functions we propose is 

WK^2^^^K{\n^\). (46) 

We denote by (A^,) the weighted asymmetries obtained by substituting the functions 
(46) into eq. (3). Taking into account eqs. (44) and once more eqs. (7), (10) and 
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(22), we get 
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(47) 




where 



^ J o?'P±P± • S^KiP) J o?Vp± • SI'l(p)/ [z, (Px - zp^f] . (48) 



The hnear system (47) can be solved with respect to Cl{x), provided we assume the 
dominance of some fragmentation mechanism [2 7], which reduces the sum over / to a 
single term. Alternatively, if data relative to asymmetry for tt^ and to X-mesons are 
simultaneously available, the system (47) can be written as 



where F runs over the final hadrons tt"*", tt" and K. This new system can be solved 
if we make some assumptions, so as to reduce the sum over / to three or less terms. 
For example, we may consider separately small and large x. In the latter x-interval 
we may neglect sea contribution and the system is overdetermined, since / runs over 
two flavors. For small x, where the sea prevails, we may solve the system by assuming 
a relation between quark and antiquark distributions, in such a way that / run over 
three flavors. Once the coefficients c{(x) have been determined for any / and for a 
sufficiently large number of L, we may determine /i{ thanks to the second eq. (44) 
and to the second eq. (21). The method we have just suggested is somewhat similar 
to the purity method used by HERMES in sphtting the single flavor contribution [46] 
in longitudinally polarized distributions (see also refs.[31, 32, 47, 48, 49, 50]). 

The weighted cross sections (A^) wash out the unwanted J J contribution, there- 
fore they are particularly suitable in an experiment, like HERMES[11], where the 
target is longitudinally polarized. 

3.3 Spurious effects 

Asymmetry (26), as well as any asymmetry connected to a SIDIS experiment, is 
sensitive to processes which have nothing to do with the distribution we want to 
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(^D E^y W^f(^) =^EE4^^l(^)5{(x), 



(49) 
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determine. Such processes, whose contributions may be perhaps reduced, but not 
ehminated by means of kinematic cuts[31], are essentially of two kinds: 

(i) interactions[30] and/or correlations[32] between the active quark and the spec- 
tator partons; 

(ii) fragmentation of the target remnant diquark or decay of a light cluster [31]. 
Gluon exchange between the active quark and the spectator partons is demanded 

by gauge invariance[34]. A gluon emitted by spectator partons may interact with the 
active quark a) before ov h) after being struck by the photon. The amplitude for such 
a process interferes with the amplitude for the active parton to interact with the sole 
photon, giving rise to a "dynamic" twist- 3 contribution [2 8, 29] (see also refs. [30, 51]). 
In our case this term must be T-even, since T-odd functions [30, 51] are involved only 
in single spin asymmetries. In particular, the typc-a^ contribution corresponds to 
5fj.(a;, p_l), whose integral over the transverse momentum is f?^!^)- Calculations [29] 
within the model proposed by Qiu and Sterman[34] assure that such contributions 
are about 10% of the parton model term. 

The spurious fragmentation effects (ii)[31] and quark-spectator correlations [32] 
are nonperturbative and can be taken into account according to a strategy similar 
to the one followed by Kotzinian[31]. The method consists of a simulation program 
based on the Lund fragmentation model. The simulation [52, 53] reproduces the data 
of unpolarized SIDIS, but exhibits deviations from the purity assumption [46, 47, 
48, 49, 50] that the final hadron observed come exclusively from the active quark 
fragmentation. 

Such deviations are suitably described by fracture functions [2 6]. In our case the 
asymmetry (26) has to be modified into 

A'(|k|;Q,i.;P) = ^ (50) 



Here M.^ and 5M.^ are respectively the unpolarized and polarized fracture functions, 
which depend on z and 11^. As in the case of longitudinal polarization[31], the 
cross section for reaction (4) can be calculated by means of the simulation program, 
assuming parametrizations (37) to (39) for the nonperturbative functions involved. 
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and using some model calculation for /if [2, 20, 44]. Since the program is endowed 
with a pointer, which shows the origin of the final hadrons[31], one can determine the 
fracture functions by selecting the spurious events that result from the simulation, 
and fitting them by means of a suitable parametrization, e. g., 6Ai-^ = C6Q-^ , C 
being a constant. 

Incidentally, in the case of longitudinal polarization, these spurious effects might 
be responsible for the discrepancy between the recent HERMES results[46] and ex- 
pectations from the so-called "spin crisis" [31]. 

It is worth noticing that the weighted asymmetries wash out in part such effects, 
for example the light cluster decays, which are isotropic. 

4 Comparison with other methods 

Extraction of hi from SIDIS single spin asymmetry is affected by drawbacks similar to 
those described in subsect. 3.3. In particular this asymmetry derives a contribution 
from the Sivers effect [33], which could be even greater than the term we are interested 
in. In fact the Sivers function is chiral even and therefore is coupled to the usual 
t.m.d. fragmentation function; this is considerably larger than the chiral odd 
and T-odd fragmentation function coupled to hiT in the single spin asymmetry (see 
sect. 1). Moreover there are suggestions that this single spin asymmetry, usually 
classified as twist 2 [25], could be instead twist 3 [54, 26], as well as asymmetry (26). 
Indeed, this latter is twist 3 (see subsect. 2.2), while it was considered twist 2 in 
previous papers[16, 23, 21]. The reason for this discrepancy is explained in appendix. 
There we discuss a factor //~^, introduced for dimensional reasons, where ^ is an 
undetermined energy scale. This is usually assumed equal to M[l, 16], whereas a 
correct normalization of the t.m.d. distributions as probability densities demands 
11 = Q/l. Lastly, in our method, the complication of inferring hi from a number of 
weighted asymmetries, as described in subsect. 3.2, is largely compensated by the 
fact that giT is coupled to the ordinary fragmentation function of the pion, whose 
parameters are well determined, contrary to the Collins function. 

Among the other methods proposed for extracting transversity, we recall proton- 
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proton collisions, with muon pairs or hadrons (pions or lambdas) in the final state[55], 
and SIDIS with a lambda in the final state[14, 15]. But asymmetries in proton-proton 
collisions have been shown to be rather insensitive to /ii[55]. Concerning SIDIS, string 
fragmentation into a lambda is much more difficult than pion fragmentation; this not 
only makes statistics poorer, but also increases the drawbacks described in subsect. 
3.3[26]. 

5 Conclusion 

To conclude, Sq± may be coupled to a chiral-even fragmentation function, in particular 
to the twist-2, unpolarized t.m.d. fragmentation function of the final hadron observed. 
This suggests an alternative, convenient method for extracting hi, circumventing 
some of the usual drawbacks that plague determination of transversity. Specifically, 
we propose to measure the double spin azimuthal asymmetry in semi-inclusive pion 
leptoproduction. For reasonable values of (4 to 10 GeV'^), and under the most 
favourable kinematic conditions, the asymmetry is estimated to be at least ~ 10%. 
The suggested experiment could be performed at facilities like CERN (COMPASS 
coll.) and DESY (HERMES coll.), where similar asymmetry measurements are being 
realized or planned. As a last comment, our results confirm the crucial role of the 
intrinsic transverse momentum[16, 17, 38, 40] and of polarized SIDIS experiments [47, 
48, 49, 50] in extracting quark distributions. 
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Appendix 

Here we establish some relationships among transverse momentum dependent 
(t.m.d.) distributions, for which we adopt the notations by Mulders and Tangerman[16] 
(MT) (see also Ralston and Soper[l] (RS) and other authors[17, 21]). In particular 
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we show that, in the chiral hmit, gir equals hir {Sq± in the present paper). First of 
all, we write the correlation matrix in the QCD parton model. Secondly we compare 
it with the general expression of the correlation matrix. As a result we find some 
relationships among t.m.d. distribution functions, which turn out to hold true even 
taking into account quark-gluon interactions and renormalization effects. We give 
also an alternative proof of the relationship between qit and hiT in the chiral limit. 
Lastly, we shortly discuss our result. 



A.l - Correlation matrix in QCD parton model 



A. 1.1 - The correlation matrix 



We define the correlation matrix as 



$(x,px;P,5) = j dp- j 



{P,S\,p{y)m\P.S), 



(A.l) 



assuming the light cone gauge = 0. Here P and S arc, respectively, the four- 
momentum and the Pauli-Lubanski (PL) four- vector of the proton, with S'^ = —1; ip 
is the quark field, p is the quark four-momentum, whose transverse component with 
respect to the proton momentum is p±, i. e., 



p = p'^n++p±+p n_. 



n+-p± = n_-p±-- 



0. 



P ^ P+n+ + p-n-, (A.2) 
n+ ■ — 1, (A. 3) 

(A.4) 



Moreover x — p'^/ P'^ is the light cone fraction of the quark momentum. 

We take a frame such that the proton has a large momentum, say V, much greater 
than its rest mass M, and it is transversely polarized. Moreover we choose the z-axis 
along the proton momentum and the |/-axis along the proton polarization, so that 



P = (VM2 + p2^0,0,P), 
n± = -^(1,0,0,±1), 

P± = (Pl,P2,0), 



-5 =(0,0, 1,0), (A.5) 
P± = (0,p±), (A.6) 
|px|=0(M). (A.7) 
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Prom now on this frame will be called P-frame. For the sake of simplicity, we exclude 
T-odd terms in the correlation matrix (A.l). 



A. 1.2 - QCD-improved parton model 

The correlation matrix for free, on-shell[23] quarks in a transversely polarized 
proton reads 

E QT{x,p^)l{^ + m){l + 2T^,^,). (A.8) 

T=±l/2 

Here m and p are, respectively, the rest mass and the four-momentum of the quark, 
such that = m^. 2TSq is the quark PL vector, with — —1. qT{x,p±) is the 
probability density of finding a quark with its spin aligned with {T — 1/2) or opposite 
to (T = -1/2) the proton spin. Now we define a quark rest frame[38] - to be named 
g-frame from now on - whose axes are parallel to those of the P-frame. In this frame 
we have Sg — S'^^ = = (0, 0, 1, 0). Decomposing S'^'> = S into a transverse and a 
longitudinal component with respect to the quark momentum, we get 

= S = E^cose' + lysine'. (A.9) 

Here 

sinO' — sin9sin(f), sinO — -^^4^, sin(f) = f — , (A. 10) 

IpI Ip-lI 

p ~ {p±,xP), iy={0,t), E^ = (0,n), (A.ll) 

t = {sin9cos4>,sin9sin4>,cos9), (A- 12) 

n = {cos9cos(t),cos9sin(f), —sinO). (A. 13) 

In order to calculate Sq in the P-frame, we perform a boost along the quark momen- 
tum. This boost leaves invariant and transforms v into p/m, where 



p={\p\,Eqt), Eq^^m? + p\ (A.14) 

As a result we get 



Sa = S + 



m 
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sin9'. (A. 15) 



Here we have defined 



m 



1 



We stress that it is essential that the relative momentum of the quark with respect to 

the proton is nonzero; otherwise x would be fixed and equal to m/M[38]. Substituting 
eq. (A. 15) into cq. (A. 8), and taking into account the definitions (A. 11) and (A. 16) 
of u and 6, we get 

+ ^5q4x,pi_)^^i^^[^J]+^sine' - A + mB^ + 0(V-'). (A.17) 
Here we have set 

Qi^, p1)^ Y1 ?T(a:,P±) 5q±{x,p±)^ ^ 2TqTix,p±), (A.18) 

r=±l/2 T=±l/2 

A = E,^[fi'^,fi'_]stne', (A.19) 



B 



+ |C (l - - + ^[r^;, fi'.]^sine'. (A.20) 
Moreover we have exploited the relations — — l/2[^,;^] — p ■ S, p ■ S — p± ■ S and 
75( 



= :^(n', -n'_). 



A. 2 - Parametrization of the correlation matrix 

Now we parametrize the correlation matrix of a transversely polarized proton 
in the T'-frame. We normalize this matrix in such a way that the leading twist 
contribution of $ for an unpolarized proton, when integrated over p^, coincide with 

the product of the usual density matrix times the distribution of an unpolarized 
quark. Taking into account the operators of the Dirac algebra[23, 16, 24], we get 

^± = ^0a + ^06 + ^l + ^2. (A.21) 



Here 
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+ 



(e + QTi,^) , 



(A.22) 



+ 



(A.23) 



$1 



(A.24) 



while $2 contains terms of 0{V We have set 



X±^ -S ■ p±/ii 



(A.25) 



and /i is an undetermined energy scale, which we shall fix below. Moreover the 
distributions involved, for which we have used the notations of MT, are functions of 
the Bjorken variable x and of the intrinsic transverse momentum p^. The term 



is interaction independent. This is evident for $oa, which consists of twist-2 operators. 
We shall show that also $ob shares this feature, although the operators involved are 
classified as " twist-3" . 

A. 3 - Relationship between hir and qit 

A. 3.1 - Comparison with the QCD parton model 

We equate the coefficients of the independent Dirac operators in eqs. (A. 17) and 
(A. 21), taking into account the first eq. (A. 2), which in the P-frame reads 



= 



$0a + $06 



(A.26) 




(A.27) 



We get 




(A.28) 
(A.29) 
(A.30) 




(1 — ei)sin9'Sq±, 
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fxX^hr = (1 - ei)sin0'EgSq±. (A.31) 
\±9iT = {l-e2)sin9'6q^, (A.32) 
X±gT = {l-e3)sin9'Sq±. (A.33) 

Here ti = m/ Eg, €2 = m/2xV and €3 = m/2|p| are the correction terms to the chiral 
Umit, which are generally small for light quarks. The terms of order O [{m? + pi)/P^] 
have been neglected. As regards ji, RS and MT have set it equal to M. We require 
the various functions to be normalized, in the chiral limit, as 5q^, which is a difference 
between two probability densities. Therefore we assume 

Ax = sine', {AM) 

which, according to eqs. (A. 10) and (A. 25), implies ^ — |p|. In particular, this choice 
leads to the relationship 

giT = (1 - e2)hiT- (A.35) 

In the chiral limit, sinO'hiT = sind' qit is the average helicity of a quark in a trans- 
versely polarized proton. 

A. 3. 2 - Equations of motion 

Owing to the Politzer theorem[37], the Fourier transform of $_|_ must fulfil the 
equation of motion (e.o.m.) for a Dirac particle interacting with the gluon field. We 
set $_L = + Since the term (eq. (A. 8)) fulfils the Dirac equation for a 
plane wave, the e.o.m. implies that depends on the quark-gluon interaction and 
is of order gV~^[l'&, 56], g being the strong interaction coupling constant. But 
includes the term $0, corresponding to eq. (A. 26) with the constraints (A.28)-(A.33). 
On the other hand, $_l — $0 includes the interaction dependent term is 
orthogonal to Tr $o(*^'± — *^'o) = 0. Therefore $0 is interaction independent and 
the relationships (A. 28) to (A.33), although deduced from the naive parton model[23], 
hold true even after inserting interactions. In particular the term $o6; although made 
up with twist-3 operators, is interaction independent and should be classified as a 
"kinematic" higher twist term. 
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The Politzer theorem survives renormahzation and off-shell effects[37], which, on 
the other hand, do not cause any mixing among the mutually orthogonal Dirac op- 
erators into which we have decomposed the correlation matrix. Therefore the rela- 
tionships (A. 28) to (A. 33) have a quite general validity; in particular they hold also 
when QCD evolution is taken into account. 



A.3.3 - An alternative proof of relationship (A. 35) 



Eqs. (A.21) and (A.l) imply 
having defined the functional 



{p, s\mi^f>'-^{y)\P, s) 



n[p{y)]^Jdp-J -^,e^''p{y)- 



Similarly, 

xVhiT 



1, 



1 



Tr{$o75[^,?i-]} = 7^ (P, ^|^(0)75[^, ?i-]^(l/)|P, ^) 



(A.36) 



(A.37) 



(A.38) 



Now we decompose the quark field into the eigenstates of a canonical representation, 
in which the quantization axis is taken along the proton polarization: 



(A.39) 



Here ip' and ip" denote respectively the "good" and "bad" component of the quark 
field, i. e., 



1 



^;y)(l/) = -(l±75^)?iH-?i-V'(y), 



(A.40) 



The field may also be decomposed into chirality eigenstates, i. e., 

1 



i'{y) = i'R{y) + i''L{y) + r{y), 

Then 

V'(O)75?i-'0(y) = 

1- 



i^'R(L){y) = 7(1 ± i5)i/i+fi-i^{y)- (a.41) 



V2 
1 



;^i^{0M^,i-My) = +^[^\Ho)my)-^[\o)^[{y) 



(A.42) 
(A.43) 
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Substituting eqs. (A. 42) and (A. 43) respectively into (A. 36) and (A. 38), we get 



(A.44) 
(A.45) 



Here we have set 



qR 



1 



n 



2V2xV 



n 



{p,s\^im'M\P:S) 

{P,S\ilj\\0)ilj\{y)\P,S) 



(A.46) 
(A.47) 



and we have defined analogously and q^. Eqs. (A.46) and (A.47) imply that 
^fl(L) probability of finding the quark in a chirality state, whereas q^^-^ is the 

probability for a quark to be in a state of the canonical represention defined above. 
In this representation the spin density matrix of a quark in a transversely polarized 
proton reads 

P = gf|T)(T|+?fU)ai, (A.48) 

I I (I) > denoting the two states of the canonical representation. On the other hand, 
the helicity operator can be written as 



A=l+)(+|-|-)(-|, 



(A.49) 



where |±) denote the two helicity states. Therefore the average helicity of a quark in 
a transversely polarized proton results in 



(At) 



rr(pA)=gf [|(+|T)r-|(-|T)P 
~ sin6' = hiTsinO' . 



i(+ii)r-i(-ii)r 



(A.50) 



But (Af|-) equals — q^ in the chiral limit. Then 



ql-ql- (?f - ql)sine'. 
The definitions (A.46) and (A.47) and equations (A.42) and (A.43) imply 



(A.51) 



{p,s\moi^{yms) 



for 



m = 0, 



(A.52) 
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having set 

= ^^\ii.-^sine'[^,iiJ\Y (A.53) 

Therefore, in the chiral hmit, relationship (A. 35) holds true independent of renor- 
mahzation and off-shell effects. 

A. 4 - Discussion 

The result wc have just found appears in contrast with the behavior of the ordinary 
distributions. Indeed, gi{x), a chiral even function, is independent of hi{x), which is 
chiral odd. Therefore an experiment for determining gi{x) automatically excludes the 
possibility of inferring hi{x) and vice-versa. But transverse momentum attenuates 
this difference. A quark polarized perpendicularly to the proton momentum has a 
nonzero helicity if its transverse momentum is different from zero. Quantitatively, in a 
transversely polarized proton, a massless quark, whose average spin component along 
the proton spin is /?.it(3:^, P±), has an average helicity sin6'hiT{x,p±), which implies 
hiT{x, p_i_) = Qirix, p±) in the chiral limit. The fact that hiT is associated to a chiral 
odd operator and Qit to a chiral even one simply means that the same function can be 
deduced from different types of experiment: either from single polarization, exploiting 
the Collins effect [25], or from double polarization, as suggested in the present paper 
and by Kotzinian and Mulders [21]. 
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